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The artificial separation of a full-theory mode into distinct collinear and soft modes in SCET leads 
to divergent integrals over rapidity, which are not present in the full theory. Rapidity divergence 
introduces an additional scale into the problem, giving rise to its own renormalization group with 
respect to this new scale. Two contradicting claims exist in the literature concerning rapidity scale 
uncertainty. One camp has shown that the results of perturbative calculations depend on the precise 
choice of rapidity scale. The other has derived an all-order factorization formula with no dependence 
on rapidity scale, by using a form of analytic regulator to regulate rapidity divergences. We deliver 
a simple resolution to this controversy by deriving an alternative form of the all-order factorization 
formula with an analytic regulator that, despite being formally rapidity scale independent, reveals 
how rapidity scale dependence arises when it is truncated at a finite order in perturbation theory. 
With our results, one can continue to take advantage of the technical ease and simplicity of the 
anal3rtic regulator approach while correctly taking into account rapidity scale dependence. As an 
application, we update our earlier study of WW production with jet-veto by including rapidity scale 
uncertainty. While the central values of the predictions are unchanged, the scale uncertainties are 
increased and consistency between the NLL and NNLL calculations are improved. 


I. INTRODUCTION 

The paradigm and framework of effective field theory 
(EFT) is essential in our conceptual understanding of 
quantum field theory as well as in practical calculations. 
In an EFT, by definition, we have removed (i.e., inte¬ 
grated out) the modes of the full theory that are necessar¬ 
ily highly off-shell with the given initial and final states, 
so that the Lagrangian only contains the modes that can 
be on-shell and may produce singularities in loop and/or 
phase-space integrals. Also fundamental in EFT is to 
have small parameters besides coupling constants and as¬ 
sign to each field in the theory a definite power-counting 
rules in terms of those parameters, so that we have mani¬ 
festly well-controlled expansions at the Lagrangian level. 
Often, defining such power-counting rules requires clas¬ 
sifying the EFT modes into smaller groups of modes and 
introduce a separate field for each group. 

Soft collinear effective theory (SCET) [TH5] [3 [7] is 
an excellent example of EFT with such further mode 
separations. In a typical application of SCET to col¬ 
lider physics, we have collinear modes, i.e., the energetic 
modes nearly parallel to one beam direction, and anti- 
collinear modes, i.e., those nearly parallel to the other 
beam. It is also often convenient or necessary to intro¬ 
duce soft mode^which are between collinear and anti- 
collinear modesQ To introduce our notation and nomen¬ 
clature, let us quickly review the power-counting rules 
for those modes. Let n and h be the 4-vectors (1,0, 0,1) 
and (1,0,0, —I), respectively. Any 4-vector v can then be 
decomposed uniquely as n^/2 + v- n^/2 + 

where = fi-v and V- = n-v, and n-v± = n-v± = 0. 


^ We restrict our discussion to SCETn-type observables, for which 
collinear, anti-collinear and soft modes suffice. Other observables 
may require different modes such as ultra-soft modes. 


Then, a collinear momentum p, an anti-collinear momen¬ 
tum p, and a soft momentum ps are characterized by the 
following scaling behavior in powers of a small power¬ 
counting parameter A ^ 1: 

P = {P+,P-,P1-) ^ (1, A)M, 

p= (p+,P_,P_l)--(A^1,A)M, (1) 

Ps = {ps+,Ps-,Ps±) ~ (A,A,A)M, 

where M is the characteristic high energy scale of the 
process in question. Perturbative computation in SCET 
is a dual expansion in A and the coupling constant Og. 
All the three modes have the same level of virtuality, 
^ p 2 ^ ^ A^M^, but their rapidities rj (along 

the beam axis) are vastly different, i.e., ^ A“^, ^ A, 

and ~ 1, respectively. As noted earlier, the EFT should 
describe those three groups of modes by three separate 
fields (j){x), (fix), and (fsix) (even if they all originally 
belong to the same field $(a;) in the full theory) so that 
we can assign a definite power of A to each field, and 
thereby each term, in the Lagrangian. 

Such separation of modes, however, can lead to ad¬ 
ditional divergences in loop and/or phase-space inte¬ 
grals [8HT0]. Just like we get ultraviolet (UV) divergences 
by ignoring the boundary (or cutoff) between the highly 
off-shell modes we have integrated out and the nearly 
on-shell modes we have kept, we get rapidity divergences 
by neglecting the boundaries between collinear, soft, and 
anti-collinear modes. Then, just like UV divergence leads 
to a renormalization group (RG) with respect to a scale 
parameter p, rapidity divergence leads to its own RG 
with respect to a different scale parameter v [TTVmj . As 
usual, the p RG equations (RGEs) can be used to resum 
large logarithms arising from a (Lorentz invariant) hier¬ 
archy of scales in the problem. Similarly, the zz-RGEs 
can be used to resum large logarithms from a (frame- 
dependent) hierarchy of rapidity scales (as in Eq. 0 ). 
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A natural question, then, is whether there are ad¬ 
ditional scale uncertainties associated with the precise 
choice of initial and final values of v in solving the v- 
RGEs, analogous to the familiar scale uncertainties from 
the choice of /i. At first sight, it might appear that the 
i^-RG should not constitute an independent source of un¬ 
certainties. Firstly, both p- and i/-RGEs resum the same 
large logarithms of the form log(l/A), because the theory 
actually contains only one large ratio of scales, 1/A, as 
can be seen in the scaling laws Q. Since there is only 
one kind of logarithms being resummed, one might argue 
that there should only be one source of scale uncertain¬ 
ties. Secondly, recall that the usual p scale dependence 
arises because Og depends on p and we do not know its ex¬ 
act p dependences That is why physical quantities like 
cross-sections depend on p even though we demanded 
that they should not. In contrast, being a Lorentz in¬ 
variant quantity, as cannot depend on v, because v is 
a rapidity cutoff, i.e., a proxy for the frame-dependent 
boundaries separating collinear, soft, and anti-collinear 
modes. Then, the i/-independence requirement for physi¬ 
cal quantities would just render the perturbation series to 
be trivially independent of order-by-order in Og, leaving 
us no additional scale uncertainties. Indeed, in Refs. m- 
[n], all-order factorization formulae were derived that are 
explicitly v independent, based on a form of analytic reg¬ 
ulator [IS] for cutting off rapidity divergences. However, 
in a seminal paper on this subject m, a thorough, ex¬ 
plicit analysis demonstrates that the i^-RG constitutes 
an independent source of scale uncertainties, where v is 
varied independently of p. A recent work |19j also shows 
that rapidity scale dependence can have significant nu¬ 
merical impact. 

In this paper, we provide a simple resolution to this 
puzzle. We will show exactly how v dependence is con¬ 
cealed in the seemingly ^-independent factorization for¬ 
mula in Ref. m, by deriving an alternative all-order fac¬ 
torization formula. Our formula is equivalent to that of 
Ref. [17] if we were to compute to all orders in pertur¬ 
bation theory. They are not equivalent when truncated 
at a finite order, however, and our formula vividly shows 
how it acquires v dependence at a finite order and allows 
us to quantify the scale uncertainty associated with the 
precise choice of u. Our understanding also captures how 
v dependence arises in spite of Og being v independent, 
as well as how it constitutes an independent source of 
scale uncertainties even though there is only one kind of 
logarithms being resummed. 

For the derivation of our factorization formula, we will 
employ a form of analytic regulator to regulate rapidity 
divergences that is similar to, but different from, what 
has used by Refs. iniiiHi. Our analytic regulator shares 
the same great advantage that there is no need to per- 


^ In practice, there is also scale dependence from parton distribu¬ 
tion functions (PDFs), but that is beside the point. One could 
talk about QED instead of QCD. 


form ‘zero-bin subtraction’ |20j to avoid double-counting 
modes in the overlapping regions, as the integrals in 
the overlapping regions become scaleless and thus van¬ 
ish once consistently multipole-expanded m- This type 
of analytic regulator is also shown to preserve fundamen¬ 
tal properties such as gauge invariance [T8|. The dif¬ 
ference between our and their forms of analytic regula¬ 
tor is that our regulator treats the collinear and anti- 
collinear sectors in a symmetric manner by introducing 
two independent rapidity RG scales, while the regulator 
of Refs. [T71[TS] treats the two sectors asymmetrically but 
needs only one rapidity scale. As a consistency check, we 
will show in Appendix that our factorization formula 
can also be obtained from the regulator of Refs. [HIIH]. 
Therefore, we can continue to take advantage of the enor¬ 
mous technical simplicity and ease of the analytic regu¬ 
lator approach for practical calculations while correctly 
taking into account rapidity scale uncertainties at the 
same time. 

As an application of our factorization formula and 
method of estimating rapidity scale uncertainty, we will 
update our earlier work of the WW production with jet- 
veto at the LHC [52], which used the formulation of 
Ref. [T7] and thereby lacked scale uncertainties from ly 
dependence as well as suffered from another accidental 
cancellation that led to an underestimation of the total 
scale uncertainties. We will see that, while the central 
values of the predictions remain completely unaltered, 
the scale uncertainties increase and the consistency be¬ 
tween the next-to-leading logarithm (NLL) and next-to- 
next-to-leading logarithm (NNLL) calculations improves. 

II. THE SET-UP, INGREDIENTS, AND GOAL 

Even though the conceptual points we will be mak¬ 
ing are quite general, for definiteness we discuss a SCET 
formalism that is applicable to the resummation of jet- 
veto logarithms for the production of a color-neutral ob¬ 
ject [23], such as a higgs boson [TT] 1231121] [H] |55H57] . 

W^W~ [221111] or other di-bosons [22], at a hadron col¬ 
lider. Jet veto rejects an event if it contains a jet whose 
transverse momentum is larger than a prescribed jet-veto 
scale Using M to denote the invariant mass of 

the color-neutral system, and keeping only the relevant 
features explicit, we can express the differential jet-veto 
cross-section with respect to M a^ 

^H{p) Zsip,v,P) B{p,iy) B{p,v), (2) 

where /x is the usual RG scale associated with UV diver¬ 
gence, while both u and u are rapidity RG scales. For 

^ We consider the simplest case where is constant, but it 

could in principle be chosen to depend on other variables in the 

process, such as the rapidity of the jet m- 
^ If the color-neutral state is a higgs boson, it would be a total 

jet-veto cross-section with M fixed to the higgs boson mass. 
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generality, we have introduced two rapidity scales be¬ 
cause in principle the location of collinear-vs-soft bound¬ 
ary is independent of that of soft-vs-anti-collinear bound¬ 
ary. We now define and explain each of H, Z^, B and B 
below. 

The hard function H (p.) is the squared amplitude for 
the hard parton collision with the center-of-momentum 
energy M that produces the color-neutral system in ques¬ 
tion. Being identical to the corresponding squared am¬ 
plitude in the full theory, the hard function has no ra¬ 
pidity divergences and is hence independent of v and v. 
The scale /i has the usual interpretation as the scale or 
boundary beyond which the modes would be regarded as 
necessarily highly off-shell. So, /i is a proxy for the level 
of virtuality of the degrees of freedom in the EFT, which 
suggests /i^ ~ ^ ~ Ps ^ A^M^. Further noting 

that p_L ~ p_L ~ Ps_L ~ AM ^ p™*° in the presence of 
jet-veto, we expect p ~ Later we will see how this 

choice is indeed forced upon us by the formalism. The job 
of p-RGEs is thus to take us from the high scale p ^ M 
(the hard scale) where the SOFT is matched onto the full 
theory, down to the low scale p ~ p™*° (the factorization 
scale) where the cross-section (§ is evaluated, resum¬ 
ming some of the logarithms of the form log(M/p™*°). 
(The remaining logarithms of the form log(M/p™*°) will 
be resummed by rapidity RGEs, as we discuss below.) 

The beam functions B{fj,,u) and B{fj.,D) [JT] |32] are 
respectively the PDFs for the collinear and anti-collinear 
partons with the jet-veto. That is, if we schematically 
write a quark PDF as 

( 3 ) 

where ip and N are respectively the collinear quark field 
and the nucleon in question, and the summation goes over 
all possible hadronic states X, then the beam function for 
the same quark and nucleon is given by 

R^5]'(iVHX)(XHiV), (4) 

v 

which is the same as before down to every little (im¬ 
plicit) detail except that the summation only goes 

over X that satisfies the jet-veto condition. The collinear 
beam function i?(p, v) suffers from rapidity divergence 
when we send fc_|_ (the -I- component of integration vari¬ 
ables) toward zero and thus ‘invade’ the territory of soft 
modes. The regulator we introduce in Section m reg¬ 
ulates this divergence and trades it for a scale v. This 
regulator does not introduce any other new scale, so we 
expect the logarithms from fc_|_ integration to have the 
form log(z^/^P+), because the only physical scale in the 
-|- direction is fP+ (where f is the parton momentum 
fraction and P+ /2 the proton momentum in the collinear 
beam). This suggests that the good scale choice that 
minimizes the logarithms should he v fP-e- Similarly, 
the anti-collinear beam function R(p, P) has rapidity di¬ 
vergence as k- —> 0, which is regulated in favor of a scale 


D. Our regulator by design treats the collinear and anti- 
collinear sectors symmetrically, so we expect P ^ 
(where f and P_ /2 are the parton momentum fraction 
and proton momentum in the anti-collinear beam). Inte¬ 
grations over the T component are also divergent, which 
we regulate by dimensional regularization acting on the 
d — 2 transverse dimensions with the scale p. As we al¬ 
ready discussed, the T component controls the degree of 
virtuality of the partons, so we expect the p-dependent 
logarithms to have the form log(p/p™‘°) in the presence 
of jet-veto, suggesting the choice p ^ 

Finally, Zs(fj.,i',P) is called the soft function, which is 
a matrix element of a SOFT operator consisting only of 
soft modes, evaluated between states containing only soft 
modes. Rapidity divergences arise when we send or 
fc_ to infinity and invade the land of collinear or anti- 
collinear modes. The regulator trades the k+ and k_ 
divergences for v and P, respectively. Alternatively, in 
our context, the expression ([^ tells us that Zg can also 
be regarded as a renormalization constant that absorbs 
the rapidity divergences and v- and t'-dependences of the 
product of operators BB. In the former interpretation, 
the soft modes are ‘integrated in’ in the SGET and we 
can directly calculate the soft function diagrammatically. 
In the latter, the soft modes are not in the theory and 
Z^ is determined by matching ZgBB onto the full theory 
at some appropriate scale. Either way, what we actually 
do in practice are essentially identical, while both inter¬ 
pretations can be useful conceptually. So we will discuss 
the two viewpoints in parallel below. Most importantly, 
whether it is viewed as a soft function or renormaliza¬ 
tion constant, Z^ should not know anything about f.P+, 
fP-, nor M. Thus, the degree of virtuality imposed by 
jet-veto, p ~ i® the only physical scale entering Zs, 

so the good scale choice must be given by i/ ~ ^ p to 

avoid large log(://p) and log(f'/p). 

We can now state our goal. For B and B, we have high 
scales ly ~ fP+ and P ~ fP-, which are both of order 
M. For Zs, we have low scales iz ^ P ^ fi, which is of 
order The job of rapidity RGEs is then to take 

us from the high scales to the low scales, resumming the 
logarithms of the form log(M/p™*°), the same form as 
what the p-RGEs resum. As we alluded in Section 
this is not surprising as there is actually only one large 
ratio of scales, 1/A, in the theory. This is also reflected 
to the correlation between the rapidity scales and the 
virtuality scale, v ^ P ^ ji. So, our problem is to analyze 
how imperfect this correlation can be at a finite order 
in perturbation theory and to quantify the associated 
uncertainty. 


III. (RAPIDITY) RG EQUATIONS 

In this section, we will derive an all-order form of fac¬ 
torization formula to be used in Eq. ([^, that is, the 
precise form of “Zs(p, t/) R(p,;/) i3(p, i/)” there. The 
reader who wishes to skip the derivation and read about 
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the discussion of rapidity scale uncertainty may directly 
move on to Section IHl 

To regulate rapidity divergences, we employ the follow¬ 
ing form of analytic regulator, where for each phase-space 
integral^ with momentum k we insert 

(^) + (^) 


out), the only dependence on besides through as{fJ,) is 
via the following logarithms: 

= log (pvet0)2 > ^-L = ^^g ^-veto)2 ' 

A. The /i-RG equations 


For a collinear k and at the leading order in A (i.e., 
0(A°)), we should neglect k- in comparison with in¬ 
side 0(fc+ — fc-), which renders the step function trivial. 
Thus, this regulator in practice simply amounts to in¬ 
serting {y/k+)°‘ in the phase-space integration for each 
collinear particle, and similarly just inserting {yjk-)'^ for 
each anti-collinear particle. This symmetric treatment of 
the collinear and anti-collinear sectors then implies that 
B and B are identical up to the trivial relabelling v 
a <->■ d, etc. 

As all regulators do, the regulator (§ introduces artifi¬ 
cial new scales into the problem, v and which are prox¬ 
ies for the boundaries separating collinear modes from 
soft modes, and anti-collinear modes from soft modes, 
respectively. The requirement that physical observables 
should be independent of v and v leads to rapidity RGBs. 
We also introduce separate jet-veto scales and 
for the collinear and anti-collinear sectors, respectively, 
as they can in principle be chosen independently and it 
has proven useful to do so Hug 

Now, as we take the a —>■ 0 and a —>■ 0 limits, l/a and 
1/d poles appear. Just like a 1/e pole in dimensional 
regularization comes with log/i, the 1/a and 1/d poles 
are respectively accompanied by 

^ ^ 

Here, ^P+l‘^ and fP_/2 are respectively the energies of 
the collinear and anti-collinear partons that enter the 
hard function iJ, while T+/2 and P_/2 are the energies 
of their respective parent hadrons. These logarithms are 
the only sources of v and v dependences. The form of 
L can be understood by observing that the only physical 
scale entering in integration in the collinear sector is 
The form of L then follows from symmetry. 

The beam functions also depend on the scale fi of 
dimensional regularization, which regulates divergence 
from k± integration. Then, since the jet-veto scale is 
the only physical scale in the T direction in the EFT (as 
the hard physics of order M has been already integrated 


^ As shown in Ref. El, rapidity divergences only occur in phase 
space integrals. 

® However, we still assume ~ fo avoid introducing an¬ 

other hierarchy of scales into the theory. We will see later that 
the formalism would indeed break down if is very different 

from 


Without loss of generality, the /^-RGEs can be written 
as 


log P(/i, tts, Lj.) LP , 

d 


p—Q 

oo 


( 8 ) 


/r— log B{p, = tts, T_l) PP , 

^ p=0 


where the /r dependences are in as, T_l, and L_l, while the 
v and V dependences are in L and L, respectively. The 
functions /1 p 1 in the B equation are the same functions 
as those in the B equation, owing to the symmetry of 
analytic regulator 

The equations are actually too general and their 
forms can be constrained significantly. Physical observ¬ 
ables of our interest, such as the WW jet-veto cross- 
section, depend only on the product Z^BB, as in Eq. ([^. 
Being a renormalization constant, cannot depend on 
dynamical variables such as Pf, or P_. Alterna¬ 
tively, being the soft function, i.e., a matrix element of 
soft fields between soft states, it cannot depend on the 
momenta of the (anti-)collinear partons nor of the collid¬ 
ing hadrons. We therefore require that the sum of the 
right-hand sides of the /r-RGEs (|^, i.e.. 


OO OO 

A = ^ /(P)(C, a„ L^) PP + ^ /'P)(e, a„ L^) PP , (g) 

p—0 P—0 


should be independent of P+, and P_ at any val- 
ues of p,, V, and v. These four variables, however, are 
not all independent of each other but are constrained as 
^P+^P- = because we are considering the quantity 
of the form ([^, which is a function of M. So, we instead 
consider 

where £ is a Lagrange multiplier varied independently 
of f, P+, and P_, so that the condition dAldi = 
0 gives us the constraint back. Then, the require¬ 
ment dA/dP+ = 0 implies 

/(')=£, = = = (11) 


The first relation here means that f^^\^,as,L±) is actu¬ 
ally independent of i.e., Lj_) = 

because the Lagrange multiplier i is by definition inde¬ 
pendent of the remaining independent variables, P+, 
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and P-. However, ^ may still depend on other param¬ 
eters and variables. Most generally, it may depend on 
Og. It cannot actually depend on L±_ nor L±_ for the fol¬ 
lowing reason. Our symmetric treatment of the collinear 
and anti-collinear sectors implies that if ^ depended on 
L_l, it should also depend on L± symmetrically, i.e., 
L±, L±) = i{as, L±, L±). However, the above con¬ 
dition f^''-\as,L±) = I tells us that £ does not depend 
on Lx, so by symmetry it does not depend on Lx ei¬ 
ther. We thus conclude that L±) = f^^\as). 

Then, applying this and the conditions ( [lT| ) to the re¬ 
quirement dA/d^ = 0, we get 


^/(°\e,ag,Lx)=0. 

We thus have narrowed down the ^-RGEs Q to 
log v) = L -h Lx), 

/i^ log H(^, v) = L + Lx). 


( 12 ) 


(13) 


Next, since the hard function H{fj,) is independent of v 
and V, the v and D dependences in the product BB should 
be completely cancelled by in the product ZsBB. So, 
we must have 


d 


logZs(Ai,z/, v) 


where 


-/lil(ag)Ls+F(o)(ag,Lx,Lx), 

(14) 


vv 

Ls = log ^ . 


(15) 


Cancelling the remaining ^ dependence in \og{ZsBBH), 
we then get 

Lj_) + /(o)(ag, Lx) + Lx, Lx) = -4j. 

(18) 


B. The i/-RG equations 


Physical quantities must be independent of the renor¬ 
malization scheme. In particular, the solution of RGBs 
must be independent of the shape of the path in the /i- 
ly-P space and can only depend on the endpoints of the 
path. This path independence requires that 


d d d d 

log log S(/r, v) 


d^i dv 


= , 


(19) 


where the second equality has used Eq. (13). Thus, 


^d/d^, acting on v dlog Bjdv gives us a independent 
function, f^^\as). So, if there are any i^-dependent terms 
in V d log B/dv, they must be independent of /x so that 
they get annihilated by fid/d^. However, since the per¬ 
turbation series is an expansion in Og, every term does 
depend on /x via Og. Therefore, there cannot be any in¬ 
dependent terms in v dlogB/dv. We thus arrive at the 
general form of rapidity RGBs: 


d 

v-^logB{ij.,v) =g(/x), 

d 

v—\ogB{ij.,v) =g(/x). 


( 20 ) 


We can determine and a combination of /1°1 and 
^(0) fj-om the requirement that the /x dependence in the 
product ZgBB must be cancelled when Z^BB is multi¬ 
plied by the hard function iJ, i.e., fid\og{ZsBBH)/dfi = 
0. The /x-RGE of the hard function can be parametrized 
as 

/^^logi?(/x) = 2rLM+ 47 , (16) 


where Lm = log(M^//x^) with the coefficients T (the 
cusp anomalous dimension) and 7 (the anomalous dimen¬ 
sion) that can be calculated order-by-order in Og using 
perturbation theory. To identify we demand that 
the derivatives and M d/dM commute when act¬ 

ing on \og{Zf,BBH). This identifies the M dependence 
in Lm with that in L -|- L from Eq. (13), i.e.|^ 


/(il(ag) = 2 r. 


(17) 


Mathematically, a /^-independent integration constant should be 
included in the right-hand side of the relation 0 . However, 
such integration constant is absent in perturbation theory as all 
terms in RGEs contain Qs, which does depend on /i. 


where 


00 00 

9{l) = LP_, g{fi) = '^g^P\as)L^j_. ( 21 ) 

p—0 P—0 


(Note that these are expansions in L± and L± unlike the 
/i-RGEs (^, which were expansions in L or L.) We have 
written g^ as a function of only Og, as it turns out that 
g^P^ is independent of ^ (or ^) as we will see shortly. We 
again require that the product Z^BB to be independent 
of IX and P, which implies that the rapidity RGEs for Zs 
must be given by 


J^^log^s(M,G v) = -g{p), 
_ d 

v-^\ogZs{g,v, ix) = -g{p). 


( 22 ) 


The path independence for Zg does not provide us 
with any new constraints. Eirst, the Zg rapidity 
RGEs (22) trivially lead to [ix (3/9ix, ix cl/clix] log Zg = 
0. Also, the path independence condition (19) for 
B and the similar condition for B readily imply that 
[jjLd/dg, vd/dv] logZg = [gd/dg,, Pd/dP] logZg = 0. 
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C. The all-order factorization formula and the 
choice of scales 


We are now ready to solve the rapidity RGBs (20) 


and (22) for an arbitrary but fixed /r. Solving them is 


trivial as the right-hand sides of all the i/-RGEs are in¬ 
dependent of V and v. We solve for B starting from an 
‘initial’ point v = to an arbitrary final v, B from 
to D, and Zs from to Then, combining all 

the solutions, we get 


Zs(/x,i/, ly) B{p,v) B{fj.,iy) 


(23) 


with 


There are two things to note here. Firstly, the right 
choice of the factorization scale at which we should eval¬ 
uate the factorization formula ( [2^ is on the order of the 
jet-veto scale p™*”, rather than the hard scale M. Sec¬ 
ondly, our formalism would not be applicable to the case 
in which and are hierarchically different. In 
such case, there would remain large logarithms of the 
form log(p™‘°/p™*’°) that are not resummed by our for¬ 
malism. 

The remaining scales Vs, Vu, and Pg should be cho¬ 
sen to minimize L, L, and Lg. These logarithms vanish 
if 

= ^P+ , = ^P- , Vs = i>s= (30) 

The factorization formula (1231 then becomes 


ZsBB{^) = Zs(/r, i/g. Vs) J^b) S(p, v^). (24) 

To determine what value of /i we should choose when 
we evaluate Zs{ij,,v,v) B{ii,v) B{fi,D) in the factoriza¬ 
tion formula ( |23[ ), let us examine more closely the path 
independence condition Applying the relation 0 
to it, we can rewrite it as a condition for g{p): 


= f^^Kas) = 2T. 


(25) 


Satisfying this condition order-by-order in L± in the ex¬ 
pansion (21), we get the following nontrivial recursion 
relations among the coefficient functions g^P\as)'. 


2g(i)+/35'(°) =/W = 2 r, 

=0 (p>2), 


(26) 


where g'^P^ = dp^^Vd^s and /3 = pdog/d/r. Let us 
parametrize g^^\ T, and /3 as 


g^^^ = tts di + Ug d2 + ■ ■ ■ , 
r = as(ro -l- OsTi -(-•••), 
(3 = —2Q:sas(/3o + Os/3i -I- ■ 


(27) 


•), 


where Og = as/(47r), and all the coefficients, di 2 ,.... 
To,!, ., /3o,i,- -) can be calculated in perturbation theory. 
Since none of these coefficients depend on the functions 
g^P'^ are indeed independent of ^ for all p as we alluded 
earlier. Combined with the parametrization (27), the re¬ 


cursion relations (26) give 


— OsTo + Q-s (Ti -I- )3odi) -)-••■ 

5 ( 2 )= +■■■ 


(28) 


Note that g^P^ = 0{aP) for all p > 1. This means that 
in order for us to be able to reliably truncate the expan¬ 
sions (21) at any finite order in Og, the logarithms L± and 
L± must be small, i.e., L±,L± <C This is possible 
only if we choose 


Zsip,v, v)B{p, v)B{g,v) 

_/ f V'Y K Wb fig, , (31) 

-\iPj (fpj 

However, since our formalism is valid only for p™*° ^ 
p™*°, there is no practical benefit to be gained by choos¬ 
ing the two jet-veto scales differently in experimental 
analyses. So, hereafter, we will mostly focus on the spe¬ 
cial case. 


PT‘o=P^?to, L^=L^, g = g. 


(32) 




, „veto 

Pt 


Pt ; 


The factorization formula (31) now takes a very simple 
form: 

/ „2 ^ ^ ^ 

Zs{p,v,v)B{fi,v)B{g,v)=\-^\ ZsBB{g). (33) 

Before closing this section, we would like to remind 
the reader that the only properties of the analytic regu¬ 
lator (^ that were actually used in the derivation of the 
formula ( [^ are the forms of the logarithms ^ and ([^ . 
Therefore, any other regulator that only produces those 
logarithms will lead to exactly the same form of factoriza¬ 
tion formula as Eq. (331. It is actually even more robust 
than that as we will see in Appendix where a regula¬ 
tor that gives different forms of logarithms still arrives 
at essentially the same factorization formula, although 
it seems difficult to pin down exactly how robust it is. 
Finally, let us also remark that it should be straightfor¬ 
ward to repeat our analysis for objects other than the jet- 
veto beam and soft functions, such as the transverse mo¬ 
mentum dependent parton distribution functions (TMD- 
PDFs) for small px resummation. 


IV. MAIN RESULTS AND DISCUSSIONS 

Let us now examine what happens to the factoriza¬ 
tion formula 0 or (|33|) if we deviate from the scale 


(29) choice (30). We are now ready to show that, while such 
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deviations cancel out order-by-order in as, the cancel¬ 
lation is incomplete if the series is truncated in a way 
consistent with the EFT power-counting rules, rendering 
rapidity RG to be an independent source of scale uncer¬ 
tainty. 


A. Finite-order truncation and additional scale 
uncertainties from rapidity RG. 

If we could calculate ZsBB exactly or to all orders 
in perturbation theory, the scale choice (30) would be 


unique and unambiguous. This is a corollary of the fact 
that the right-hand sides of j^-RGEs (20) and (22 1 are 
independent of v and D. For example, one might ob¬ 


ject to the choice (30) that the vanishing of only im¬ 
plies VsVs = so we could instead choose Vs = rii 
and Vs = r“^/i with some r. Glancing back at Eq. (23), 


we see that this (apparently) alternative choice would 
amount to multiplying the prefactor 
in the factorization formula (31) by a factor of while 
changing Zg to Zs(fj,,rfj,,r~^fi) at the same time. But 
the j^-RGE (22) tells us that Zs(/r,r/i,= r~^^^Zs, 
so r cancels out. In fact, any one of t'e, 7 'b, and Vs 
can be varied independently, e.g., Vj^ = iP+ —>■ r^P+, 
without producing any net effect. This ‘r invariance’ is 
exact, and since Og is independent of v and it is exact 
order-by-order in Og. 

However, the whole point of using an EFT is to re¬ 
sum the large logarithm log(l/A). So, in the EFT power 
counting, aglog(l/A) is by definition parametrically re¬ 
garded as 0{1). Thus, the EFT perturbation series is 
not a literal power series in Og but an effective power 
series in Og with log(l/A) counted as 0{a~^). When the 
factorization formula (33) is truncated at a finite order 


according to this EFT power counting, the r invariance 
is no longer exact. Let us say we are aiming at an 0{a^) 
accuracy with the EFT power counting (i.e., the so-called 
N^+^LL accuracy). Then, the right-hand side of factor¬ 
ization formula (33) should be truncated as 


2 - 

^^1 ZsBBU) 

L J computed to 0 ( 0 ^) 

. 

where hn are Og-independent coefficients defined via the 
expansion g{^) = a". Recalling the defini¬ 


tion (24) and the scale choices (29) and (30), we see that 


the object [ZsBB{iP)'\ has no logarithms oi v or v at all, 
nor any large logarithms of /i. That is why the object 
[ZsBB{fi)] is calculated above to literally 0{a^), with 
no need to take into account the EFT power counting. 
In contrast, the exponent of is computed above 

to C>(as+^), because \og{M^/g?) ~ log(M^/(p™*°)^) ^ 
log(l/A) ^ 0{a~^) and the EFT power counting must 
be applied. 

Let us now examine how Z^BB changes from the ‘cen¬ 
tral value’ (34) if we deviate from the scale choice (30). 


Again, looking back at Eq. (23) (with g = g) and us¬ 
ing the i/-RGEs (20) and (22), we see that the pref¬ 
actor (^ 2 ^^ 2 ^/ 1103 - 1 -/ 1205 -I-'+?!.fc+ia^+^ gg|. multiplied 
j,/ti 03 +/t 20 ^H h/»fe+io(+^_ Pqj, product of ma¬ 
trix elements \ZsBB{g?)\^^ there are two possi¬ 

bilities depending on how the r dependence is calcu¬ 
lated, which agree to 0(a^) but differ at and 

hence leading to different estimates of rapidity scale un¬ 
certainty of One way is to start with an ©(ob¬ 

truncated “initial” Z^BB evaluated at the scale (30) 


and then evolve it to a different scale choice by us¬ 
ing the ©(Og )-truncated //-RGEs. The other way is 

to directly evaluate Z^BB at ©(a^) at a point away 
from the choice ([30|). In the first method, the prod- 


to 0(a") 

. Then, combined 


uct of matrix element \ZsBB[g)\^^ would be sim¬ 
ply multiplied by - 

with the factor of coming from 

(/iVM 2 )^i“=+^ 2 “=+-+^'=-n«?+\ we would be left with a 
simple net r dependence, . This is not optimal 

for the purpose of estimating rapidity scale uncertainty, 
however, because it is only sensitive to a single coefficient 
hfc+i, which may be accidentally small or large. On the 
other hand, the second method leads to 


Zs{g,,v, v) B{p,v) B{g.,v) 

/ 2 - 

= {^\ 


mV 


/tict3-|-/i2Ct^H-|-/ifc+ia3+^ 


(35) 


X r 


^—h\Ots — hoCtt - huOit 


ZsBB{g)] 


to O(a0 


which is more robust as it involves not only all oi hi 
through hk+i but also the whole perturbative series of 

ZsBB{g) to 0{a^). Hence, we expect that varying r 
by an ©(1) factor in the expression (35) should give 


us a plausible estimate of rapidity scale unc erta inty of 
©(aV^)- We therefore propose the formula (js^as our 
N^+^LL truncation of the factorization formula (33) with 


rapidity scale uncertainty appropriately taken into ac¬ 
count. 


B. Comparison with the literature 

We can now discuss exactly what is missing in the all¬ 
order factorization formula of Ref. which exhibits no 
rapidity scale dependence. They used an analytic regula¬ 
tor different from Eq. ([^, which treats the collinear and 
anti-collinear sectors differently, but as we will show in 
Appendix]^ it nonetheless leads to a factorization for¬ 
mula essentially identical to Eq. (33): 


Zs{g)B{g,,v)B{g.,v) =[^\ ZsBB{fi) 


(36) 


where ZsBB{ii) = Zs{g,) B{g,,v^) B{g,,v^) and we have 
taken as we did for Eq. (33). Then, repeat- 





























ing the analysis of Section IV A we arrive at a correct 
N^+^LL expression 




„h,. 


fc + l“a 


2 x/iiQbH- 




ZsBB{^) 


which is essentially the same as Eq. (35). 


J to 0(a’’) 
(37) 


To obtain the form of factorization formula of Ref. [H] 
from Eq. we first set r = 1, since Ref. has 

made the scale choice i^b = iP+ and Db = 

(see Appendix 0 implicitly and did not consider vari¬ 
ations from it. We also rewrite as p? IM"^ = 


[(p™‘°) 2/M2] e^^ 
would get 


Then, instead of the formula (37), we 


Zs{p)B{p,v)B{p,v) 


M 2 


Z'sB'B'ip) 


to 0{a^ 


where Z'^B'B'{p) = eS^^ZsBB{p). 


(38) 

This is the form 


presented in Ref. [IT]. Since there is no r dependence, 
there is no rapidity scale uncertainty to talk about. 
One might argue that the lack of the r-dependent pref¬ 
actor r^k+ia^*^ could be (retrospectively) justified by 
saying that it is 0{a^^^) for r = 0{1). However, 

the factor r^k+ia^'^ 


IV A 


in our analysis of Section 
is a deviation in the prefactor from the central value 
(/i^/M^)^'^“=+^^“s“' hft.fc+ias''' ^ which is an 0{q’^) quan¬ 
tity in the EFT power counting. This “fundamental am¬ 
biguity in the exponentiation of the rapidity logarithms” 
was first pointed out in Ref. m- This ambiguity itself 
should be regarded as part of rapidity scale uncertainty, 
and our form ([37| (or (35)) captures this ambiguity. 


V. AN APPLICATION: THE WW 
PRODUCTION WITH JET-VETO 


Let us now apply the factorization formula (35) to the 
WW production with jet-veto at A: = 1, i.e., at NNLL. 
For the exponent coefficients hi and / 12 , we need 0{a^) 
calculations. Combining the parametrization (27) and 


solutions (28) gives 


hi — y—(di + ToLj.), 
47r 


/l 2 = 


(47r)2 


<^2 + (hi + Podi)L± + 


ro/3o 


(39) 


For the WW production at the 7- and 8-TeV LHC runs, 
the qq initial states dominate, so the coefficients di 2 and 
Top must be found for the fundamental representation 
of SU(3). From explicit diagrammatic calculations, one 
finds di = 0, while ^2 can be found in Ref. [2T], where 
it is denoted by while Top and (3o can be found in 

Ref. [13]. The dependence of jet-veto on the jet-radius 


parameter R enters the calculation through and 

this R dependence was first studied in Refs. [33ll35] [2T] . 

For [ZbBB{p)'^, we need 0{as) calculations. From an 
explicit 1-loop calculation, one finds that 13 is given by 

/■i 

-8=51 / (40) 

Jf ^ 
i^q,9 ^ 


where (/>^ is the PDF of parton i for the proton and 

<5(1 - z) Sj 


L 2 

1 - as( +7oAp 






(41) 


where and TT-^p-i can be found in Appendix A of 

Ref. |28|. B is given by the same formula except for the 
obvious replacement of “g” by “g”. Then, the whole 


ZqBB can be determined by matching it onto the corre¬ 
sponding full QCD expression at p = and evolving 
it to another p by using the /r-RGEs (13) and (14). Alter¬ 
natively, Zs can be directly calculated diagrammatically 
by regarding it as the soft function. Either way, one finds 


Zb = 1. 

We can now update the scale uncertainty estimates for 
the WW jet-veto cross-sections presented in Ref. [22] . 
Note that the central values of the predictions, defined 
by r = 1 and p = are clearly unaffected. How¬ 

ever, scale uncertainties will necessarily be modified. In 
particular, the scale uncertainties arising from varying 
r from 1/2 to 2 will be added in quadrature to those 
from varying p from p5j?*°/2 to 2p™*°, treating the two 
sources of uncertainties as independent. We will not in¬ 
clude power corrections as they were shown to be very 
small in Ref. [36] for the WW production at the 7- or 
8 -TeV LHC run with = 25-30 GeV. For other pro¬ 
cesses and/or choices of parameters, power corrections 
can become important, and a formalism for smoothly 
transitioning from the resummation-dominated regime to 
the fixed-order regime is described in Ref. m- Including 
the uncertainties from non-perturbative effects [35] [25] 
is beyond the scope of this paper. 

In Fig. [^ the NLL and NNLL WW jet-veto cross- 
sections are shown for the ^/s = 8 TeV LHC run with 
the anti-fcx jet algorithm |39j with the jet-radius param¬ 


eter R = 1. For the NLL, /12 was discarded and B and 
B are set to the corresponding PDFs. The resumma¬ 
tion’ [4(11142] is included in all the cases as in our earlier 
work [22] The plot on the right is based on the naive 
factorization formula as used in Ref. [22], where we see 


® For a thorough discussion of scale uncertainty associated with 
‘vr^ vs no vr^’, see Ref. m- 
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FIG. 1. The NLL and NNLL jet-veto cross-sections for WW 
production for the 8-TeV LHC run with R = 1. The scale- 
uncertainty bands in the left plot are obtained from the cor¬ 
rect factorization formula ( |35|), while in the right-hand plot, 
they are obtained following |22l. 


note that a modest change in /12 can lead to substantial 
change in the final numerical results as /12 appears in the 
exponent of the ratio /r/M. Not having /i 2 , on the other 
hand, the NLL bands of the i? = 1 and i? = 0.4 cases 
have a similar width. 


ACKNOWLEDGMENTS 


PJ thanks Patrick Meade for discussions on scale un¬ 
certainties in the WW jet-veto cross-sections. We also 
thank our anonymous referee for pointing out to us that 
our previous version of Eq. (35) was only sensitive to the 


single coefficient hk+i and thus overly simplistic as dis¬ 
cussed in Section |IV A TO is supported by the US De¬ 
partment of Energy under grant DE-FG02-13ER41942. 



FIG. 2. Same as Fig. but for R — 0.4. 


that the NNLL uncertainty is ridiculously small (< 1%), 
much smaller than the well-established NLO scale uncer¬ 
tainty of the inclusive WW cross-section, ^ 3-4% [43]. 
On the other hand, the plot on the left is obtained from 
the correct factorization formula (35), which now exhibits 
scale uncertainties of reasonable size. 

In Fig. we have repeated the computations for i? = 
0.4, a typical value chosen in LHC experiments. Clearly, 
the plot on the left, which uses the correct factorization 
formula (35), displays a better convergence of the NLL 
and NNLL predictions than the plot on the right, which 
is the one presented in Ref. f^ . 

Comparing the two plots on the left, we see that the 
NNLL band in the i? = 0.4 case is much broader than 
that in the i? = 1 case. This broadening is mainly due 
to a larger value of /12 in Eq. (35). The large value of ft .2 
for i? = 0.4 comes from a large ^2 due to a logi? multi¬ 
plied by a large coefficient enhanced by factors of Ca and 
TpUf (further multiplied by numerical factors of 0(10)) 
arising from gluon splittingj^ Perhaps more importantly. 


® For a detailed study of It dependence at Ofaf) and the resum- 
mation of log/?, see Refs. m and m, respectively. Including 
the results of those studies into our analysis is beyond the scope 
of the present paper. 


Appendix A: An analytic regulator with a 
collinear-anti-collinear asymmetry 


Instead of the analytic regulator ([^, we can also use a 
different analytic regulator where we use 



for both the collinear and anti-collinear sectors iniiiH], 
which can be quite convenient for practical calculations 
due to its simplicity. While the asymmetric treatment of 
the collinear and anti-collinear sectors means we cannot 
obtain B from B by simply relabelling as —)■ 

etc., this asymmetric regulator has the simplicity that it 
only introduces one rapidity scal e, u. Since the action 
of the asymmetric regulator in the collinear sector 
is the same as the symmetric case ([^, the form of L is 
unmodified from the form (§. On the other hand, in the 
anti-collinear sector, the presence of {v/k+)°‘ instead of 
(i>/fc_)“ implies that L should be changed to 


L' = log 



(A2) 


This form can be understood from the fact that k+ = 
k^/k- (with fcx = |^_l|), where the fcx integration is cut 
off by dimensional regularization at the scale fi while the 
fc_ integration picks up ^P-, the only physical scale in 
the anti-collinear sector. 

Then, repeating a similar analysis as we did with the 
symmetric regulator in Section |III[ we find that the fi- 
RCEs (13) and (141 are modified as 


log B{fi, v) = /l^l(as) L -t- A_l) , 

logi3(/i, ly) = -/l^l(as) Z' -b Zx), (A3) 

k-^^ogZsikw) = F^°\as,Lj_,Lj_) , 
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where the B equation is the same as before as nothing 
has changed in the collinear sector, while the B and 
equations are changed. Due to the lack of symmetry, the 
function is not a priori the same as and Bg has 
no reason to be the same as either. Also observe 
that the right-hand side of the Zg equation above has no 
dependence on The constraints 0 and ( |18[ ) remain 
as before with the obvious replacements —>• 

/o(as,l^) and B(o)^Fm. 

For rapidity RGBs, we find that the z/-RGEs (20) 
and are modified as 


d 


log S(/r, v) = ^ g^^\ag) L\ 


_L ’ 


p—Q 


log s(/r, v) = -^ g^^Kas) LP_ , 


(A4) 


p—0 

oo 


’^^^ogZsifJ., I') = - '^g^p\as) {LP_ - ). 

p=0 

Again, the B equation is unmodified. The negative sign 
in the B equation here is a co nsequence of the negative 
sign in the B equation in (A31 combined with a relation 


in the g-v space. 


similar to Eq. (25) from path-independence requirement 


Solving Eq. (A4) for B from Vg to an arbitrary final 


i^, B from i/e to v, and Zg from Vg to ly we find that the 
factorization formula (|23l) is modified as 


(A5) 


Zg{n,v)B{ii,v)B{ix,y) 


with ZgBB{fj,) = Zg[iJL,Vg) B{^,v^) B{^,v^). We can 
again show ‘r invariance’ as we did in Section [IV A[ and 
the minimization of tells us that the choice for in 


Eq. (30) should be modified to 


iP- 


(A6) 


while z/b remains unchanged. For /i and Vg^ we again see 

Vg g 


by matching ZgBB to the correspond¬ 


ing full QCD quantity, since ZgBB (by definition) only 
depends on /i, i^g, and while the QCD counterpart 

has only fi and Alternatively, if we integrate in 

the soft modes in the SCET and regard Zg as the soft 
function, it is clear that is the only physical scale in 
the problem and hence jj, ^ Vg ^ Hence, we choose 


Vg = g, as in the symmetric regulator counterpart (30). 


For the case of practical interest we have 

g = g, which makes Vg cancel out in Eq. (A5), and the ra¬ 


tio t'e/r'e becomes as in the symmetric regu lato r 

case. Furthermore, from the iz-RGE for Zg in Eq. (A4|, 


we see that Zg becomes v independent when g = g (while 
it may still depend on /r). Therefore, the factorization 


formula (A5) reduces to Eq. (36), which is essentially 
identical to the symmetric regulator counterpart (0 . 


The discussion of rapidity scale uncertainty in Section IV 
then goes through essentially unmodified and leads to the 
final factorization formula (37). 
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